Support of Non-separable Multivariate Scaling Function. 

Irina Maximenko * 



C^T ' Abstract 

o 

We make an estimation of the support of a multivariable scaling function for an arbi- 
trary dilation matrix. We give a method of calculating the values of the scaling function 
^ . on a tight set using the knowledge of the size of the support. 

Q 

m 

^ ! 1 Introduction 

■ For construction of wavelet bases it is often convenient to use a scaling function, a function 

which satisfies the following functional equation 

tp( x ) = 2^2c q ip(2x - q), xeR. (1) 

The equation (1) is called a scaling equation and the sequence of coefficients {c q } is called 
a mask. The scaling function can be built from the mask. Finite support masks are of 
the most interest because they generate wavelet functions with compact support. The 
knowledge of the support, particularly, allows to apply an algorithm of constructing the 
scaling function on a tight set. In one-dimensional case if function (p satisfies equation 

OO ! N2 

<p(x) = 2 c q f{2x -q), x G R, (2) 

q=Nl 



O 



then the support of scaling function is contained between Nl and N2 (you can see, for 
example, [2], [1]). 

In d-dimensional case the coefficient in scaling equation is a matrix which satisfies some 
natural demands. Dilation matrix is an d x d integer matrix such that all its eigenvalues 
moduluses are more than unit. We define the norm of a matrix as the operator norm from 
R d to R d : 



\\M\\ = max \Z\Xmm*\, 
where M* is conjugative to M matrix. For dilation matrices we have 

lim j|M~"||=0 and lim \M n x\ = 0, Vx G R d , x ^ 0. (3) 

n— *+oo n^+oo 

Hence the set {M~ J k}j £ %^ k£L d x& tight in R d . Note that the norm of matrix M _1 isn't 
obliged to be less than unit. For example matrix 
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is a dilation matrix, its eigenvalues are ~ 2.3028 and —1.3028, nevertheless the norm of 
inverse matrix approximately equals to « 1.1233, that means that ||^4 _1 || > 1. 

Let m := |detM| (detM is the determinant of matrix M), then m is integer, more 
than unit. The equation 

<p(x) = m ^2 c q (p(Mx - q), x G R d (4) 

q€Z d 

is called a scaling equation and the function <p is called (M, c) scaling function. 
The cascade operator T c is the linear operator given by 

(T c f)(x):=mY,c q f(Mx-q). (5) 

q£Z d 

The iteration scheme T£f = T C (T™ -1 /), n = 1,2, ... is called a cascade algorithm. As a 
initial function Fq we can take any peacewise continuous function with compact support 
which satisfies the following conditions 

^ Fq(x + q) = 1, F(0) = 1 and F(u) is continious in the origin. (6) 

q£Z d 

For example as a Fq we can take the characteristic function Fq : = X[-i/2,i/2) d or = 
d 

]7 B(xj) where B(xj) = 1 — if \xj\ < 1 and B(xj) = in another case. 
k=i 

Let's define 1-periodic by each variable function 

mo(u) := c q e- 2ni ( q ' u \ u G R d . 

Obviously if {c q } has finite support then function mo is a trigonometric polynomial. The 
scaling equation in frequency domain is 

(p( u ) = m (M*- 1 u)(p(M*- 1 u). 

Repeating this procedure (suppose that Fourier transform of scaling function is continuous 
in origin and ^(0) = 1) we get 

oo 

(p(u) = Y[m (M*- j u). (7) 

i=i 

For convergence of the infinite product it is necessary that mo(0) = 1 or (the equivalent 
condition) 

J2 c k = 1. (8) 

If mo is a trigonometric polynomial then the infinite product converges uniformly on 
compact sets. 

Let the space S' be a set of tempered distributions defined on a set of test functions 
S = S(R d ) (infinitely differentiate functions which with all their derivatives decrease on 
infinity faster than arbitrary power function) . The sequence of distributions (or functions) 
f n converge in S' to distribution / if lim (f n , g) = (/, g) for all functions g G S. 

n—>oo 

Distribution / equals to zero in neighborhood U of the point x$ if for any test function 
g which is not zero only on neighborhood U we have (/, g) = (see for example [3]). 
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If distribution / isn't equal to zero in any neighborhood U of the point xq then x$ is 
called essential for functional /. Totality of all essential points is called the support 
of distribution /. The support of distribution / corresponding to usual continuous or 
piecewise continuous function / is the closure of the set on which f(x) / 0. 

In book [4] it is proved that the condition (8) is sufficient for convergence of the cascade 
algorithm in S'. So 

Theorem 1 [4]. Let {c q } be a finite mask with (8). Then scaling equation (4) has unique, 
up to constant factor, decision ip G S' with compact support. This decision is given by (7). 
More over, for any distribution with compact support f G S' sequence f n = T n f converges 
in S' to function c • if where factor c equals to /(0). 

In section 2 the support of a multivariate scaling function is estimated for an arbi- 
trary dilation matrix and for masks which provide the convergence in S' of the cascade 
algorithm. Then obviously the estimation is true for other types of convergences. 

In theory if we know the mask we can find the scaling function and the wavelet function 
because we can make inverse Fourier transform in (7). However we can use a simpler on 
practice method of construction (p. In section 3 we show how the knowledge of size of 
support can be used for calculating values of a scaling function on a tight set (in one 
dimensional case this method is described, for example, in [1]). 

2 The Estimation of Scaling Function Support 

Lemma 2 Let F±, ...,F n , ... be peacewise continuous functions with compact supports on 
M. d . Let the sequence of functions F\, F n , ... converge in S' . Suppose that supp F± C 
Ai, supp F2 C A2, ... ,supp F n C A n , where the sets A\, A2, A n , ... G R d are 
closed full spheres with the same center. Let the sequence of spheres A n converge which 
means that there exists a sphere A such that its radius is the limit of radiuses of spheres 
A n . Then the support of limit function F is in A. 

Proof. We are going to show that any point out of A doesn't belong to the support of 
distribution F. Let's fix an arbitrary point ao 6 \ A. Since the set A is closed there 
exists neighborhood U of point A which is separated from the set A. Since lim A n = A 
there exists a number N such that all spheres A n , n> N will be separated from U. Let 
test function g G S be not zero only in the neighborhood U. Then (F n , g) = for all 
n > N and hence (F, g) = 0. That means that point ao doesn't belong to the support of 
distribution F. 

Note. The lemma is true if instead of full spheres we consider rectangle parallelepipeds 
with faces parallel to coordinate hyperplanes with the same center. The convergence of 
rectangles means here the convergence of edges. 

Denote the set of indexes q such that c q is not zero. The mask c q has a finite support 
hence £1 is a finite set. Remind the definition of cascade operator (as an initial function 
Fq we take a piecewise continuous function with compact support which satisfies (6) ) 

F n (x) = (T c F n -i)(x) := mJ2c q F n -i{Mx - q). (9) 

qen 

First we estimate the support of limit (scaling) function when ||M _1 || < 1. 
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Theorem 3 Let M be a dilation matrix and ||M _1 || < 1. Suppose that mask {c q } has 
finite support and (8). Then for the support of limit function ip the following estimation 
is true 

supp ip C {x G R d : \x\ < 1 _ | M _i | }» 

where Q := max Id. 

Proof. Let's fix an initial function Fq which is a piecewise continuous function with 
compact support that satisfies (6). Denote its support as Qq. The support of function 
Fq is a compact set hence there exists a full sphere with radius R which contains Slo, 
R := max{|x|, x G ^o}- Applying cascade operator (9) to function Fo we get function F\. 
Function F\ is not zero if Mx — q G for any g 6 fior what is the same x G M~ 1 (tto+q). 
Denote the support of function F\ as Then for x G f2i it is true that 

\x\ < \M-\n + q)\ < WM^WR + \\M- l \\Q. (10) 

Denote the support of function F n as fl n . Let's prove by induction that function F n equals 
to zero outside of the set of x where 

\x\ < \\M- 1 \\ n R + QQIM^W + IIM" 1 ]] 2 + - + ll M_1 ir)- (11) 

For n = 1 the assertion is true (10). Suppose that estimation (11) is true for n and we 
prove it for n + 1. Applying cascade operator (9) to function F n we see that function 
F n+ i is not zero if Mx — q G O n for any q or what is the same x G M~ l {£l n + q). Then 
for x G fin+i it is true that 

\x\ < \\M~ 1 \\Q+ ||M- 1 ||(||M- 1 || n .R + Q(||M- 1 || + WM^f + ... + WM' 1 ]^)) = 

= WM^W^R + QiWM^W + HM" 1 !! 2 + - + HM" 1 !!^ 1 ). 

The induction is completed. By theorem 1 cascade algorithm converges in S'. Because 
||M _1 1| < 1 the limit of the right side (11) exists and is finite. Using lemma 2 and directing 
n to infinity we have 

Q IIM" 1 !] 

i-||M-i|| ' * Gsu pp^- . 

Note. The support of scaling function ip doesn't depend on the size of support of the 
initial function Fq. 

Corollary 4 Let d = 1. Then dilation matrix equals to m, m G Z, \m\ > 1. Suppose 
that mask {c q } has finite support and (8). Then 

Q , 

supp ip C {x G K : \x\ < - — : -}. 



Corollary 5 Let M be a diagonal dilation matrix d x d. Suppose that mask {c q } has 
finite support and (8). Then 

supp if C {x = (xi, ...,x d ) G R d : \x k \ < ,.9 7 , k=l,...,d}. 

- 1 
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Let now the norm of matrix M 1 be not necessary less than unit. Consider first the 
case if matrix M is an Jordan box size s with eigenvalue A, |A| > 1. 



Mx 



/AO 
1 A 
1 



V o o 



\ 





( si \ 



X2 
X3 



1 A J \x s J 



( A.r, X 

X\ + Ax 2 
x 2 + Ax 3 

V + Xx s J 



(12) 



Lemma 6 Let dilation matrix M be a Jordan box size s with eigenvalue A. Suppose that 
mask {c q } has finite support and (8). Then for each coordinate of the support of the limit 
function (p we have 



if |A| ^ 2, then \xk\ < 
*/ |A| = 2, 



Q 



l 



l 



|A| - 2 (|A|-l) fe 
then \xk\ < Qk, k = l,...,s 



k — 1, s; 



Proof. Matrix M is a dilation matrix hence |A| > 1. Let's fix an initial function Fq a 
piecewise continuous function with compact support which satisfies (6). Its support is 
contained in the sphere of radius R, R := max{|x|, x € ^o}- Applying cascade operator 
(9) to function Fq we get function F\. Function F\ is not zero if Mx — q G fio for any 
q e or what is the same x G M~ 1 (VLq + g). As above denote as VL\ the support of 
function F\. Then \Mx\ < R+\q\ < R + Q for x G fii. First coordinate of vector Mx 
equals to Axi hence |Axi| < \Mx\ < R + Q. When we multiply matrix M by vector the 
first coordinate doesn't depend on others. Then by the same reasons as in theorem 3 we 
have 

l4 n) |<^ + Q(^ + ^ I + ... + ^ r ) =:^ 1; n=l,2,3,.... (13) 
Using lemma 2 and directing n to infinity we have 

Let's prove by induction that for support of function F\ for each coordinate we have an 
estimation 

\x^\<(Q + R)^ + ^ + ... + 2^j=:A lk , k = l,...,s. (15) 
For k = 1 the assertion follows from (13) 

\x?\<(Q + R)^. 

Suppose that estimation (15) is true for k (k = 1, s — 1) and we will prove it for k + 1. 

+ Xx kli\ < \ Mx \ <Q + R- 
Using the induction postulate we have 

l4+il£^j«? + - R + ' 4 ifc) = 
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= (Q + R)[ ^7 + ^2 +- + 



,|A| ' |AP - |A|*+i 
The induction is completed. So the estimation for function F\ support is proved. We get 

(n 



function F n from the formula (9) hence if \x^ \ < A n k then 



ki-l + ^l <Q + A n -1 k 

Hence 

l4 n) l < -^(Q + A n _ lk + A nk _ 1 )=:A nk , k = 2,...,s, n = 1,2,3,... (16) 
We get a recurrent formula 

A n k = pq- (Q + A n _i k + A n fc _i) . 
Directing n to infinity we get the estimation for support of limit function if 

A-oo k = y^j {Q + ^oo k + Aqq k~i) j k = 2, s. 
From this equation for k we have 

Let's prove by induction by k that 

A ~ k = Q {\^i + Jw^ + --- + W^)' fc = 1 '-' a - (18) 

For = 1 the assertion is true (14). Suppose now that estimation (18) is true for k (k = 
1, s — 1), and prove it for k + 1. 

Ax> fc+1 = |^| _ j- (Q + ^oo fc) = 

Q + Q[ TTT—T + ,,m 1 ,N9 + - + 



|A|-i V "VIAI-1 (|A| - 1)^ (|A|-1) 
„f 1 1 1 

= Q TTl— T + TTTi TTo +- + 



|A|-1 (|A| - l)^ (|A|-l)*+\ 
this completes the induction. Calculating the sum of finite geometrical progression for 
|A| 7^ 2 we have 

for |A| = 2 we have 

4x>k = Qfc, fc = l,...,s. B (20) 
Note. By proving lemma 6 we didn't use the fact that q are integers. 

An arbitrary non-degenerate matrix M can be represented as M = C~ 1 GC ', where G 
is a Jordan matrix (consists of Jordan boxes), matrix C is unitary. Each of Jordan boxes 
(including boxes size lxl) generates a subspace invariant to multiplying by matrix M. 
Let's denote the rectangular parallelepiped P <G M. d as follows: 

if in row p of matrix G there is a simple eigenvalue X p , then \x p \ < i ^9_ 1 ; 

if in row p of matrix G there is a beginning of a Jordan box of size s which corresponds 

to eigenvalue X p / 2 then |x p+fc | < (* ~ (p^Fip ) ' fc = lj -' s; 

if in row p of matrix G there is a beginning of a Jordan box of size s which corresponds 
to eigenvalue X p = 2 then |x p +fc| < Qk, k = 1, s. 
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Theorem 7 Let M be a dilation matrix and all its eigenvalues X±, Xd be real. Suppose 
the mask {c q } has finite support and (8). Then the support of limit function p is contained 
in CP, where matrix C and set P are defined above. 

Proof. The matrix M has real eigenvalues hence in formula M = CGC~ l matrices G and 
C are real. Applying cascade operator (9) to function F n we get 

F n+1 (x) =m^2c q F n (CGC- 1 x-q), x E M d . (21) 

qen 

Making the change of variables x = Ct we get 

F n+1 (Ct) =mY,c q F n (CGt-q), t£R d . (22) 
qen 

Denoting F*(t) = F n {Ct) the formula (22) will be: 

Fl +1 (t)=mY^c q F l n {Gt-C- l q), t E R d . (23) 

Directing n to infinity we have 

Vl (t) =m^2c q ip 1 (Gt-C~ 1 q), teR d , (24) 
qen 

where tpi(t) = ip(Ct). 

Matrix G consists of Jordan boxes. Each of them corresponds to a subspace invariant 
to multiplying by matrix G. Hence using the note after lemma 6, we can apply lemma 6 
to function ipi, to each Jordan box separately. If the corresponding eigenvalue is simple 
then we apply corollary 4. Then we get rectangular parallelepiped which we denote as P. 
Then we apply matrix C to the bounds of support of function p>\ and we get the set CP 
which will contain support of function p. 

3 Values of Scaling Function on a Tight Set 

Let coefficients c q be not zero on the set £1, where Q is the finite system integer d- 
dimension vectors. If we know the bounds of support of scaling function <p, then using 
scaling equation (4), we can find values of function (p on tight set {M _:? fc} J> o i kei d - ^ we 
know the values of scaling function tp in integer points then equations 

ip(M~ 1 x) = m c q^( x ~ a )i 
qen 

<p(M~ 2 x) = m^c g ^(M _1 x - q), 
qen 



uniquely define values of function p in points Af j G Z, k G Z d . 

For finding values of function p in integer points we will use scaling equation (4) once 
more. Suppose values <p(k±), p(k^) are not zero. Making transformation of summation 
index in the right side of (4) we get 

p(k) =m^2 c q p(Mk -q) =m^2 c M k- P <p{p)- 
qen peni 
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Denote r = (ip(ki), <^(A;jv)) and denote B as the operator generated by matrix cmu- 
Then the matrix of operator B looks like 



B = m 



I c Mk!-ki CMk!-k 2 
CMk 2 -k! CMk 2 -k 2 



c Mki-k N ^ 
c Mk 2 -k N 



\ c Mk N —ki c Mk N —k 2 



c Mk N -k N / 



In the matrix form we have r = Br. The values of function (p in integer points are the 
coordinates of eigenvector r which corresponds to eigenvalue 1 of matrix B. The vector r 
is normalized so that function ip satisfies partion of unit: 



If cascade algorithm strong converges then the eigenvector corresponding to eigenvalue 1 
of matrix B is unique up to constant. 
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ip(ki) + tp{k 2 ) + ... + y{k N ) = 1. 
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